Europhysics Letters 



PREPRINT 



Dynamic Modeling of the Electric Transportation Net- 
work 

Alessandro Scire(*), Idan Tuval, Victor M. Eguiluz 

Instituto Mediterrdneo de Estudios Avanzados IMEDEA (CSIC-UIB), E07122 Palma 
de Mallorca, Spain 



PACS. 89. 75. -k - Complex systems. 

PACS. 89.75.Fb - Structures and organization in complex systems. 
PACS. 05.65.+b - Self-organized systems. 



Abstract. - We introduce a model for the dynamic self-organization of the electric grid. The 
model is characterized by a conserved magnitude, energy^ that can travel following the links of 
the network to satisfy nodes' load. The load fluctuates in time causing local overloads that drive 
the dynamic evolution of the network topology. Our model displays a transition from a fully 
connected network to a configuration with a non-trivial topology and where global failures are 
suppressed. The most efficient topology is characterized by an exponential degree distribution, 
in agreement with the topology of the real electric grid. The model intrinsically presents self- 
induced break-down events, which can be thought as representative of real black-outs. 



Introduction. - The electric grid is a critical infrastructure for our economy and society. 
Recent events, ranging from the large-scale blackouts a few years ago to the California crisis 
today [1], highlight the need to enhance the insight on the electric grid, complementing the 
traditional technological analysis [2] with new transversal points of view. Our approach is 
to study the electric grid macroscopic behavior rather than to dissect individual events. At 
that macroscopic scale, the electric grid exhibits behaviors typical of complex systems. For 
instance, on the basis of 15-years time series of transmission system black-outs of the U.S. 
electric network [3,4], it has been proposed that the electric grid may be a self-organized 
critical system, operating at or near a critical point. 

A signature of the electric grid is that it can be represented as a complex network, where 
nodes are the generators and the links the transmission lines. Recent research in complex 
networks has shown that a detailed knowledge of the topology of a communication or trans- 
portation network is essential for the understanding of cascading failures [5-8] . While some of 
these studies have focused on the topological robustness of the underlying network to random 
failures and targeted attacks [9], other research have considered dynamic processes on static 
networks [10,11]. However, these studies do not consider the network as a dynamic entity 
whose evolution is driven by the action of the nodes [12-17]. In this Paper, we present a 
dynamic model aiming to describe the growth and evolution of a transportation network. The 
network growth relies on the need of resource distribution in a heterogenous environment. 
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Fig. 1 - (Left panel) The average failure events per node (circles) and the integer part of e/A (stars) 
versus fr. Smax = 100, /i = 1/2, Nstep = 2500, N=2500. Inset. Thick line: £ = A; diamonds: ^isolated; 
solid lines: other iso-energy lines s = 5, 10, 15, respectively. (Right panel) Temporal evolution of the 
total number of failures per node, fr — 0.3. 



The Model. - We consider N dynamic elements located at the nodes of a two-dimensional 
square lattice. Each element i is characterized by its size 5^, drawn from a probability dis- 
tribution p{s). To each element i are associated two dynamic variables: the load (energy 
consumption) l\^ and the supply (available energy) We assume that (i) the load of element 
i is described by l\ = rrii^ ^f^idi where is a constant value and represents a fluctuation 
term; and (ii) initially at t = the available energy and the constant load at each element are 
proportional to its size Si 

f? = frS^ (1) 

rui rrirSi , (2) 

where fr and are constant values and in general fr > rrir. Starting from an initial condition 
where each element is isolated, the network grows as follows. At each time step: 

1. If the load overcomes the supply at element i (// < /|), a failure occurs. 

After a failure, the failing element i chooses a target neighbor through a wiring strat- 
egy. Following the empirical results observed in communication networks with spatial 
constraints [20] , the target node j is chosen such that it maximizes the function 

AiJ) = -^ , (3) 

where dij represents the Euclidean distance between the two elements, and the exponents 
/3 and 7 indicate the preference for size and proximity, respectively, in the wiring. 

2. The supply at element i is updated according to fi. 

where kj is the number of links possessed by the element j, j running over the neighbors 
of the element z, V{i). The prescription given by Eq. allows the energy to be dis- 
tributed through the links of the network as it is needed, depending on the instantaneous 
load of each linked element, making next overloads more unlikely. 
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Fig. 2 - (a-c) Network structure; d) Degree distribution corresponding to the network displayed in a) 
(diamonds), b) (circles), and c) (stars). Lines correspond to fitted curves. 



Each N steps time steps, the links are actually set and the network is consequently updated. 
The choice of the time step for the network construction is a further degree of freedom of our 
model. In real electric grids, the time scales of the fluctuating demand are much faster than 
the time scales at which the network is modified. Therefore, we choose to update the network 
each Nsteps ^ 1 time steps of the local dynamics. 

The initial total supply and the total load in the system are given by 

E° = J2frSi = NfrAf (5) 

i 

i 

where Af and are two coefficients that depend only on the geographic distribution of city 
sizes P{s). Under the transformation given by Eq. (Q, the total energy is conserved. Thus 
in the remainder we will remove the superscript. 

To model the actual geographic distribution of city sizes, the distribution of Si follows a 
Zipf law [18] of the form P{s) = As~'^ in the interval (l,5max) and A = Smax/(<§max — 1) is 
the normalization coefficient. Fluctuations are described either by the logistic map in the 
chaotic regime, or by white noise [19]. 

The average energy e and average load A per node in the system are now given by 

^ = T~r = fr r~ = /r0(<5max) (7) 
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Fig. 3 - The normalized standard deviation cr/{k) of the degree distribution versus 7 for [3—1. 
Insets: degree distribution corresponding to 7 = and 7 = 3, with fixed (5 — 1. 



A = — = mr6(Smax) + -^/rrh-^= ^^ rnr^{s^^^) + - V^*^(5max) . (8) 

iV 4 ^ -Smax ~r J- ^ 

It is worth noting that at any given time the maximum possible load in the network is Amax = 
mr6(smax) + ^/^^(smax)- Thus for a system composed by isolated nodes following the same 
dynamics, the total energy that guarantees that all the nodes have access to the energy they 
need is 

^^isolated = Amax = ^rOl^max) + V^*^(^max) • (9) 

Dynamics of failures. - The transportation network evolves as a consequence of the 
overload. As more links are added, sets of nodes will become connected forming clusters 
sharing their energy. We introduce the time average of the number of failures per node as 
an order parameter of the system. The order parameter shows a transition that depends 
on the available energy in the system, For low values of /r, e < A, the system evolves 
towards a fully connected network in which the total number of failures is of the order of the 
system size (Fig. ^). The total supply in the system is not able to sustain the needs of the 
nodes. If e > A, a statistically failure free network is formed, accounting for a cooperative 
behavior in which the network is able to effectively redistribute the energy throughout the 
system. Therefore the condition e = \ represents a transition point. We remark that, since 
the transitions occurs at a value of £ lower than ^isolated (see figCJi, inset), in our model a 
dynamic sharing of the resources in a heterogeneous environment is energetically favorable. 
However, isolated break-down events (a huge amount of elements fails at once, see Fig. 
are still present in the stable regime. These events represent global cascade failures induced by 
local fluctuation, which are amplified by propagation through the network. Global cascades 
are not suppressed, as it is not possible to satisfy the total energy needs. However, as one 
furtherly enters the stable regime, global cascades are less likely to occur, representing rare 
events that do not prevent the average number of failure to approximate to zero. 

We find interesting to analyze the transition depending on m^. If and fr are inde- 
pendent, i.e., the load and the initially available energy at each element are independent, the 
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Fig. 4 - Degree distributions for three values of fr = 0.3, 0.4, 0.7. Other parameter values in (a-c) 
as in Fig. 13 d) Cumulative failure histograms varying fr across the transition point. 



transition e = X is simply given by fr = mr-\- |y^rrv^(5max)/0(5max)- However, we should 
expect that the available energy and the load are related: elements with a high load will also 
have large amount of available energy. In this case depends on fr accounting for the fact 
that the energy availability in each node is meant to supply the local need. If, for simplicity, 
we assume that they are proportional = /i/r, with /i G [0, 1] the transition is given by 

2 (i-/i) ^('^max)/€)(-gmax) • the remainder we will assume that /i = 1/2. 



£ — 



The dynamic network. - A key dynamic ingredient is given by the dynamics of the 
network. The evolution of the dynamic network is driven by the failures of the nodes and the 
attachment function given by Eq. In the model, failures are induced by local fluctuations 
in the load, afterward propagated by the transportation network. During the transient regime, 
several failures occur simultaneously leading to the initial development of the network. After 
the transient, the network settles to a configuration where the average number of failures 
drops to zero. We characterize the emerging topologies in the stationary configurations. Two 
factors determine the network topology: the wiring strategy and the available energy e. A first 
characterization of the network topology is given by the degree distribution and the average 
shortest path length. By varying the parameters P and 7 in the wiring strategy given by 
Eq. we obtain different classes of networks. In Fig.[2lwe show the degree distribution for 
three well known topologies exhibited by our system for a fixed value of e. The fitted curves 
individuate a Gaussian degree distribution when the selection of the target element is merely 
based on its relative Euclidean distance to the failing element {P = and 7 = 1); a power-law 
link distribution when the selection of the target element j is based on its size aj {P = 1 
and 7 = 0); and finally an exponential degree distribution when the selection of the target 
element j is based both on its size aj and relative distance to the failing element {P = I and 
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Fig. 5 - (a) Average shortest path length, (b) Network efficiency. Parameter values as in Fig. 12 



7 = 3). To analyze the parameter range in which different network topologies are observed, we 
have measured the normalized standard deviation a/ (k) of the degree distribution for different 
parameter values 7 G [0,4], for P = 1. An exponential degree distribution is characterized 
by cr/{k) = 1, while it increases for a power law degree distribution. Figure El shows that 
exponential degree distributions are obtained for 7 > 2. 

Increasing the energy above its transition value = A), the degree distribution remains 
qualitatively unchanged. However, a higher number of links is required to stabilize the network 
as e approaches the transition value (Fig. This increment in the total number of links is 
due to the higher number of failures that occur during the formation of the network closer to 
the transition, as is shown by the cumulative failure histograms in Fig. 

A further topological characterization is given by the average shortest path length / [21] 

where Sij is the shortest path between the elements i and j. 

During the dynamic growth of the network, / decreases from infinity to a stable finite value. 
For a fixed E value, the minimum / is obtained for the network characterized by a power-law 
degree distribution (Fig. EK)- However, this does not mean that the network possessing a 
power-law degree distribution is the most efficient transporting the energy across the system. 
Indeed, a possible measure of the efficiency in the energy transport is the following: 

Since in our system the largest elements act as energy reservoirs, the definition of the efficiency 
77 states that the shorter is the path to the larger elements, the more efficient is the network. 
Fig. 03 shows that the network possessing an exponential degree distribution is the most 
efficient in our model. 

Discussion. - In summary we have introduced a model for the dynamic self-organization 
of a transportation network, in terms of its growth and evolution. Our model contains several 
features which are peculiar of the electric grid. The electric grid is a dynamic network that 
evolves based on local needs, failures of local elements, and energy sharing criteria. If the 
total energy is below the average consumption, the system reaches a fully connected network; 
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whereas if the system has enough energy, the emerging network reflects the cooperative be- 
havior in which the resources are globahy shared. When the cooperative behavior is achieved, 
the network stops growing and reaches a configuration with a non-trivial topology shaped by 
the wiring strategy and the available total energy. Nevertheless, local fluctuations are still 
present once the network is formed and keep affecting its functionality. Our model intrin- 
sically presents self-induced break-down events, which can be thought as representative of 
real black-outs. Another remarkable feature of our model is the heterogeneity in the initial 
energy distribution, modeling the geographic inhomogeneity in the energy availability. In a 
heterogeneous environment, the dynamic sharing of the resources is energetically favorable. 
Heterogeneity also induces the selection of a preferred topology as the most efficient in the 
energy transport, because it displays the shortest path to the bigger elements. The most 
efficient topology is characterized by an exponential degree distribution, in agreement with 
the topology of the real electric grid [22]. 
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